Achieving the Neel state in an optical lattice 
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We theoretically study the possibility of reaching the antiferromagnetic phase of the Hubbard 
model by starting from a normal gas of trapped fermionic atoms and adiabatically ramping up an 
optical lattice. Requirements on the initial temperature and the number of atoms are determined 
for a three dimensional square lattice by evaluating the Neel state entropy, taking into account fluc- 
tuations around the mean-field solution. We find that these fiuctuations place important limitations 
on adiabatically reaching the Neel state. 
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I. INTRODUCTION. 

An optical lattice is a regular periodic potential for 
neutral cold atoms [1| which enables the controlled ex- 
perimental exploration of paradigmatic ideas and mod- 
els from condensed-matter physics. This is because cold 
atomic gases generally allow for a great deal of experi- 
mental tunability. For example, Feshbach scattering res- 
onances allow for the interaction strength to be experi- 
mentally varied over a considerable range [3, l3| . Other 
quantities that may be altered include temperature, den- 
sity, and strength and shape of the trapping potential. 
In particular, an optical-lattice potential plays the role of 
the ion-lattice potential encountered in electronic solid- 
state physics. The energy bands resulting from this pe- 
riodic potential lead to a quenching of the kinetic energy 
of the atoms with respect to their interaction energy, en- 
abling the exploration of strongly-correlated phases that 
play a significant role in condensed-matter physics. 

An important model that can be studied experimen- 
tally with cold atoms is the single-band Hubbard model, 
which consists of interacting fermions in the tight-binding 
approximation. The Hubbard Hamiltonian is realized by 
cold atoms in an optical lattice when the potential is 
strong enough so that only the lowest-energy band is pop- 
ulated Q. For bosonic atoms one then commonly refers 
to this model as the Bose-Hubbard model. The theoret- 
ically predicted Mott-insulator-to-superfluid phase tran- 
sition [5| for this model has indeed been observed exper- 
imentally Q. 

The fermionic Hubbard model, referred to simply as 
the Hubbard model, is important in the context of high- 
temperature superconductivity [7|, |8| and has also been 
realized with cold atoms Q . At half filling, corresponding 
to one particle per lattice site, the ground state of this 
model is antiferromagnetic, i.e., a Neel-ordered state, for 
strong enough on-site interactions. As the filling factor is 
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reduced by doping, the system is conjectured to undergo 
a quantum phase transition to a d-wave superconducting 
state [lOl • A theoretical proof of the existence of d-wave 
superconductivity in the Hubbard model is still lacking 
and would be a major step towards understanding the 
superconducting state of the cuprates. With the recent 
experimental advances in the field of ultracold atoms, an 
experimental exploration of this issue is within reach. 

In view of this motivation, a significant problem is 
determining how the Neel state of the Hubbard model 
can be reached experimentally. In this paper, we study 
theoretically the process of adiabatically turning on the 
optical lattice [ll|, [3 ■> with the goal of determining the 
conditions required for an initially trapped balanced two- 
component Fermi gas with repulsive interactions to reach 
the Neel state in the lattice. Experimentally, the presence 
of antiferromagnetic order in this cold-atom experiment 
can be subsequently detected from shot-noise correlations 
in the density distribution [ij, [iJI ■ 

Our results are summarized in Fig. [1] For initial tem- 
peratures lower than Tp, the Fermi temperature in the 
trap, the entropy per particle in the trap depends lin- 
early on temperature as is shown by the dashed line. 
The optical lattice is then turned on adiabatically and to 
determine the final temperature of the gas we need the 
entropy per atom in the lattice. For a sufficiently smooth 
trapping potential such that the tunneling does not be- 
come site-dependent, the only effect of the trap is to place 
a restriction on the total number of particles which we 
discuss later and, other than this, we may neglect the 
trap for calculations in the lattice. Since we consider 
balanced gases here, we will at sufficiently low tempera- 
tures first enter the Mott phase with one particle per site, 
and the subsequent evolution of the gas is then described 
by the Heisenberg model for the spins alone. The result 
from the usual mean-field theory is shown for a lattice 
depth of 6.5i?R (where E-n is the recoil energy) by the 
black curve, and is equal to /cb ln(2) everywhere above 
the critical temperature Tc. Since entropy is conserved 
in adiabatic processes, the final temperature is simply 
the temperature at which the final entropy in the lattice 
equals the initial entropy in the trap. Two such processes 



arc shown by the dotted hnes for different initial temper- 
atures demonstrating that the gas is sometimes heated 
and not cooled by the lattice. Nevertheless, mean-field 
theory leads to the intuitive result that as long as the en- 
tropy per particle in the initial state is less than k^ ln(2), 
which is the maximum entropy of the Heisenberg model, 
the Nccl state is always reached by adiabatically turning 
on the optical lattice. 

The inclusion of fluctuations leads however to a more 
restrictive condition. To probe the effect of fluctuations, 
we present an improved mean-field theory which pro- 
duces a temperature-dependent entropy above Tc, as seen 
from the inset of Fig. [2l Although this approach is exact 
at high temperatures, it fails to account for spin waves 
present at low temperatures and for critical phenomena 
near Tc. By further extending the improved mean-field 
theory to reproduce the correct critical and low temper- 
ature behavior due to fluctuations, we are able to deter- 
mine the entropy in the lattice for all temperatures (red 
curve in Fig. [T|). In particular, we find that fluctuations 
lower the entropy of the atoms in the square lattice at T^ 
as 



S{T = Tc) ~ NkB ln(2) - 



3iVj2 



32kBT^{3i' - 1) 



(1) 



where v is the critical exponent of the correlation length 
f. For the case of three dimensions, z/ = 0.63 [i3|. As a 
result, the initial temperature required to reach the Nccl 
state is more than 20% lower than that found from the 
usual mean-field theory, but fortunately remains experi- 
mentally accessible. For example, with ^°K atoms and a 
final lattice depth of SE'r the Neel state is achieved when 
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FIG. 1: (Color online). The entropy per particle in the 
harmonic trapping potential only (dashed line) , in a lattice of 
depth Vb = 6.5-Er (E-r. is the recoil energy) from single-site 
mean-field theory (solid curve) and with fluctuations (dash- 
dotted curve), where Tp is the Fermi temperature in the trap. 
The horizontal dotted lines illustrate cooling and heating into 
the Neel state at constant entropy by starting in the harmonic 
trap and adiabatically turning on the lattice. 



the final temperature in the lattice is 0.0I2rF, which can 
be obtained with an initial temperature of O.GSOTf- 



II. SINGLE-SITE MEAN-FIELD THEORY. 



The Hamiltonian for the Hubbard model is given by 



n 



•^ (3 3') 



'^i.aC/ 



f^E4 



4i^J^i^Jd 



(2) 



in terms of fermionic creation and annihilation operators, 
denoted by cj ^ and Cj^a, respectively, where a labels the 
two hyperfine spin states |t) or ||) of the atoms. In the 
first term of this expression, the sum over lattice sites 
labeled by indices j and j' is over nearest neighbors only 
and proportional to the hopping amplitude given by 






(3) 



Here, Vb > is the depth of the optical lattice potential 
defined by 

V^(x) = l/o[cos2(27ra:/A) -h cos2(27ry/A) -h cos2(27rz/A)], 

(4) 
where A is the wavelength of the lattice lasers. The sec- 
ond term in the Hamiltonian corresponds to an on-site 
interaction of the strength given in the harmonic approx- 
imation by 



U = 47ra 



%Yr 



3\ 4 



(5) 



where a is the s-wave scattering length which is equal to 
174ao for '^'^K. It is well-known [15[ that at half filling 
and in the limit that U ^ t the ground state of the 
Hubbard model is antiferromagnetic and that, for k^T 'ti 
U, its low-lying excitations are described by the effective 
Heisenberg Hamiltonian 



J 



i/ = l^S,.S, 



(6) 



{3k) 



with S being one half times the vector of Pauli matrices. 
The exchange constant J = At'^/U arises from the su- 
perexchange mechanism. That is, the system can lower 
its energy by virtual nearest-neighbor hops only when 
there is antiferromagnetic ordering. 

Within the usual mean-field analysis of the effective 
Hamiltonian in Eq. ^ , the total entropy for N atoms in 
the optical lattice is given by 



5* = - 



dT 



where Fl is the Landau free energy. 



Fi^{n)^N{'^^-kBT\n 



2 cosh 



zJ|n| 
kBT 



(7) 



, (8) 



in terms of the staggered, or Neel, order parameter n = 
{—y (Sj) for the phase transition to the antiferromagnetic 
state. In the expression for the free energy, 2; = 6 is 
the number of nearest neighbors for a three-dimensional 
simple square lattice on which we focus here, k^T is the 
thermal energy, and (n) is the equilibrium value of the 
order parameter determined from 



aFL(n) 



dn 







(9) 



n=(n) 



It is nonzero below a critical temperature ksTc ~ Jz/A ~ 
(3/2) J. After solving Eq. ^ the entropy is determined 
using Eq. ([7]). The results for S and (n) obtained in this 
way are plotted as solid black curves in Figs. [T] and [2l 

The entropy S'fg of the initial normal state before 
ramping up the optical lattice is the entropy of a trapped 
ideal Fermi gas. It is most conveniently determined from 
the grand potential 



Q{^i,T) = -kBT / dep(e)ln 



,--(e-M)/fcBT 



(10) 



where /i is the chemical potential, and the effect of the 
harmonic trapping potential with the effectively isotropic 
frequency uj = [uxUjyijJzY'^ is incorporated via the den- 
sity of states p(e) = e^ /{hui)^ of the atoms. The entropy 
at fixed total particle number A^(/i) = —dil/diJ. is then 
given by S'fg = — 917/9T|^j(jv)- At temperatures much 
lower than the Fermi temperature in the trap, given by 
Tf = {3Ny^^hu;/kB, we find in this manner that ^ 
5fg = NkBTr'^T/Tp. Now, by equating the final and ini- 
tial entropies we calculate the temperature of the Heisen- 
berg spin system that results after adiabatically turning 
on the optical lattice, in terms of the initial temperature 
of the trapped Fermi gas. 

From the expression for the free energy, Eq. ([5|) we im- 
mediately see that S = Nks ln(2) for all temperatures 
T > Tc, as was shown in Fig.[TJ Although this is the cor- 
rect high-temperature limit of the entropy, temperature 
dependence will lower the entropy at Tc and therefore 
lower the initial temperature required to achieve the Neel 
state. To obtain the temperature dependence above Tc, 
we must thus go beyond single-site mean-field theory to 
include fluctuations. The simplest such model described 
below incorporates the interaction of a given site with 
one of its neighbors exactly and treats interactions with 
the rest of the neighbors within mean-field theory. 



III. TWO-SITE MEAN-FIELD THEORY. 

The two-site Hamiltonian for neighboring sites labeled 
"1" and "2" is given by 

i7 = JSi-S2 + J(z-l)|n|(S^-S^) + J(z-l)|n|2, (11) 

where the last term is a correction to avoid double count- 
ing of mean-field effects. Diagonalizing this Hamiltonian 



we obtain the free energy 



Fi^in)^N\-{z-l)J\n\' 



2/3 



■In 



2e 



-13 J/ A 



PJ 



2e'3"'/4cosh ^Vl + 4(2:-l)2|nP 



.(12) 
and find the entropy from Eq. ([7]) with the condition 
Eq. ([9|) as in the single-site model. The results are plotted 
in Fig. [21 where we see that fluctuations lower the critical 
temperature and also bring about a 2% depletion of the 
order parameter which is now less than 0.5 near T = 0. 
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FIG. 2: (Color online). The staggered magnetization (n) of 
the single (solid curves) and two-site (dashed curves) mean- 
field theories, the latter of which shows depletion at T = 
and a lowering in Tc- The entropy of both theories is plotted 
in the inset. Above Tc we see that the entropy of the two-site 
theory is temperature dependent. 



The two-site result carries the exact 1/T^ dependence 
of the entropy of the Heisenberg model at high temper- 
atures. Near T = 0, however, the entropy is still ex- 
ponentially suppressed reflecting the energy cost of flip- 
ping a spin. This exponential suppression is an artefact 
of the mean-field approximation that ignores the Gold- 
stone modes which are present in the symmetry-broken 
phase. Furthermore, critical behavior cannot be prop- 
erly accounted for by a one-, two- or higher-site model 
since, near the onset of Neel order, critical fluctuations 
extend throughout the entire lattice so one would in prin- 
ciple need to include all sites exactly. To overcome these 
shortcomings, we extend our two-site model below to all 
temperatures. 



IV. FLUCTUATIONS. 

The two-site mean-field theory produces the correct 
normal-state entropy behavior in the high-temperature 



limit, 



S{T » Te) = NkB 



ln(2) - 



3J2 



64fc2T2 



(13) 



In the low-temperature regime, the entropy is determined 
from spin-wave fluctuations prevalent near T = which 
give a black-body-like entropy, 



S{T « TJ 



iVfc, 



47r2 / ksT 



45 V2V3J(n) 



(14) 



The continuous interpolation between these two 
regimes has the additional constraint that, near Tc, we 
should obtain the correct critical behavior of the anti- 
ferromagnet, namely, the correct universal ratio of the 
amplitudes above and below the phase transition A^ /A^ 
and correct critical exponent dv — I where 



S{T^Ta) = SiTc)±A^\t\ 



diy-l 



t = (T^Tc)/Tc -^ 0=" 



(15) 
This follows from the fact that the singular part of the 
free energy density behaves as F^ /Sf", where the cor- 
relation length diverges like ^ ~ 1^1 "^ as i ^ 0. Ex- 
plicit expressions for the entropy embodying the cor- 
rect behavior in the low-, high- and critical tempera- 
ture regimes arc presented in the Appendix of this pa- 
per, and plotted as the red curve in Fig. [1] for d = 3 
using A-^ /A- ~ 0.54 and v = 0.63 [i3|, and the Neel 
temperature of T^ = 0.957fcB/J [11. Their value at T^ 
leads to the central result of this paper, namely, Eq. ([T|) 
which specifies the initial entropy required to reach the 
Neel state. 



DISCUSSION AND CONCLUSIONS. 



contribution of the magnons to the entropy diverges when 
(n) — > near Tc- One way to potentially resolve this is- 
sue is to start from the Hubbard Hamiltonian Eq. ^ 
directly but such an analysis is involved [19| and has yet 
to be carried out. 

In the above, we have focused on the d = 3 case. 
Whilst our results can easily be extended to the d — 2 
case, a more pertinent way to reach the two-dimensional 
antiferromagnet most relevant to high-temperature su- 
perconductors, would be to adiabatically prepare a three- 
dimensional Neel state, as explained in this paper, and 
then decrease the tunneling in one direction by changing 
the intensity of one of the lattice lasers. In this way, the 
three-dimensional system is changed into a stack of pan- 
cakes of atoms in the two-dimensional Neel state. Fur- 
thermore, studying doped optical lattices made by intro- 
ducing a small imbalance in the initial state may shed 
some light on the physics of high-temperature supercon- 
ductors and would be an exciting direction for future 
research. 
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As briefly mentioned earlier, there is a limit on the to- 
tal number of atoms in the trap, beyond which at low 
temperatures it is energetically more favorable to doubly 
occupy sites in the center of the trap, thereby destroying 
the antifcrromagnctic state, rather than singly occupy- 
ing outlying sites where the trap potential is larger than 
U. Thus, insisting that the system end up in the Mott- 
insulator state with one particle per site entails the upper 
bound, N < A^^ax = (47r/3)(8C//mcj2A2)3/2^ where m is 
the mass of the atoms and A is the wavelength of the lat- 
tice lasers. For ''°K atoms in a lattice with a wavelength 
A = 755 nm and depth 8i?R,, and with a harmonic trap 
frequency a;/27r = 50 Hz, A'max — 2 x 10^ which is well 
above the typical number of atoms in experiments. 

We have also attempted to determine the effect that 
fluctuations have on the entropy in a more microscopic 
manner by studying gaussian fluctuations about the 
mean-field (n) for the single-site mean-field theory in 
the low-temperature regime. But such a random-phase 
approximation has severe complications related to the 
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fact that (n) enters in the magnon dispersion as uj^ 
(|n|)|k|. Hence, as can be already seen from Eq. (fTi|) , the 



APPENDIX: ENTROPY FORMULAS 



For temperatures above Tc we use 



SjT > Te) 
NkB 



ai 



T-Tc 
T 



~T~ 



ln(2) 



with ai = 3 JV32k(k - l)klT^ and k = 3i/ - 1 ~ 0.89 
[rn | . The first term embodies the correct critical behavior 
whereas the remaining terms are present to recover the 
correct high-temperature limit. Below Tc, however, we 
have 



SjT < Tc) 
Nkv 



= -OL2 



' /Tc-Ty 




[\ Tc J 

AC(K-I)r2- 

2 T2 



where 



6 



(22 



fAn^klTJ^ 



(k-1)(k-2)(k-3) Vl35V3J3 

-ai(K-l) + /?i-ln(2)) ; 






/3i =ln2- J- 



,6{A+/A- + l) + KiK-5) 47r2fc3T3 



64Kfc|T2A+/A- 



135V3J3 



The first and last terms in S{T < Tc) embody the 
critical phenomena and allow for the continuous inter- 



polation with S{T > Tc) respectively, whereas the re- 
maining terms are included to retrieve the correct low- 
temperature behavior. 
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